Golubev and Zaikin reply: In a CommentS Aleiner, Altshuler and Gershenzon (AAG) stated that wJ "found 
that "zero-point fluctuations of electrons" contribute to the dephasing rate 1 /t v " . We made no such claims. Obviously, 
the eigenmodes of any quantum system preserve their coherence. On the other hand, single-particle properties of an 
interacting system - which are not directly expressed in terms of eigenmodes - can display a finite decoherence length 
L v in the low-temperature limit. An example is provided by a particle coupled to a harmonic oscillator battel, which 
we discuss explicitly below. This-ia also the case for the weak localization correction to the conductivity in the presence 
of interaction, where we found BI3 L v to remain finite at T ~ > 0. AAG furthermore concluded that our procedure 
is "simply wrong" because "some contributions were lost" during "uncontrollable" semiclassical averaging. This is 
not the case, as demonstrated by the fact that our approach!] fully reproduces the results of AAGB if analyzed on a 
perturbative level. The advantage of our formulation is that we can proceed beyond. „ 

The conductivity and the electron density matrix p can be expressed in terms of a path integrafl 

S ' fvrfvp exp(i£ - iS' - iS R - Si), (1) 

where 5*0 and S' Q represent the electron action on the two parts of the Keldysh contour, while iSji + Si accounts for 
J> . the interaction. We reproduce AAG's results on a perturbative level if we expand the path integral ([!]) up to first 
order in iSr + Si. To show the equivalence we perform (see Appendix A for details) some exact transformations 
starting from the formula for the conductivity of AAGB, making use of the expression for the Green-Keldysh function 
@ri r 2 ( e ) = I d r> [G> r ,{e) — G A r /(e)] [l — 2/v,r 2 ] • The result (of either approach) consists of the two groups of terms: 
G R G A G A G A and G*G A G A {1 -2p)G A . Terms of the type G R G A G R G A (claimed to be lost! in our calculation) vanish 
due to the causality principle. Thus no diagrams ocjpaths are missing in our analysis. 
t— I , The crucial difference between AAG'sB and ourotj procedures is that we do not expand in iS^ + Si, but rather 
evaluate the full path integral ([j]). The expressions Sr and Si are non-local in time and are not small at low T even 
if the interaction is weak. Hence an expansion of the path integral (jl]), which is equivalent to a Golden-rule-type 
perturbation theory in lowest order in the scattering processes, becomes insufficient at low T. The only approximation 
we employ is an RPA expansion of the effective action. This yields tractable expressions for Sr and Si, which include 
I . (within RPA)r-processes of all orders. Averaging over disorder plays no important role at this point: our resultsa and 
QO ' those of AAGH are different already before averaging. 

ON ■ The difference between AAG's and our approaches can be illustrated by the example of a quantum particle (with 
mass m and coordinate q) interacting with a Caldeira-Leggett (CL) bath of oscillators (with coupling strength 7 and 
high-frequency cutoff ui c ). Its reduced density matrix p{q\,q2) is also determined by a path integral of the type (Q). 
An expansion in the interaction iSr + Si in the long-time limit yields Fermi's Golden rule with conservation of the 
bare energy. Within this approximation, a particle with initial energy E — > at T = cannot exchange energy and, 
hence, phase coherence is preserved. This is the qualitative argument of Ref. |^, which refers to a situation where 
the total energy of the system is fixed to be the sum of energies of the nonintepacting particles. On the other hand, 
. if one does not expand but evaluates the path integral ffl) exactly, one obtainsu in the long-time limit, independent 
of the initial conditions, p{q\,q2) oc exp[— m(E)(qi — (72) ]■ Here, the expectation value of the kinetic energy of the 
• i-H ! interacting particle^t T = is (E) — 7ln(oj c /7) > 0. These results can also be derived by an exact diagonalization 
of the initial modeB; however, they cannot be ob tained from the Golden rule approach. The off-diagonal elements of 
the density matrix decay on a length scale ~ 1/ymp). Due to interaction this length scale is finite even at T = 0; it 
would diverge if one would assume (E) ~ T. The latter applies to the (obviously coherent) eigenmodes, but physical 
quantities which are expressed in the basis q will be sensitive to the decay of the density matrix p(q\, 92) as a function 
of <7x ~ Q2- We encountered a similar situation while describing the interacting electron system. A more detailed 
discussion is given in Appendix B. 

Concerning experiments we can say that our results are in a quantitative agreement with various experiments, 
especially in the quasi- 1-dimcnsional case. The disagreement reported by AAG is observed in strongly disordered 
systems with short mean free paths down to several or even one (!) Angstrom. Such systems are well beyond the 
applicability of our quasiclassical theory. An extended discussion of the experiments is given in Appendix C. 
We are grateful to our numerous colleagues and friends for support. 
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APPENDIX A: 



Here we will compare and analyze the expressions for the weak localization correction to the conductivity in the 
presence of interaction obtained in our paperp and in Ref. |[ The authors] raised an extensive critique of our results 
which includes qualitative arguments, "highlighting errors" in our calculation and comparison with experiment. The 
existence of an explicit error in our calculation would be the most serious argument of AAG against our theory. 
Therefore our first task will be to demonstrate that the AAG statement about the mistakes in our calculation is not 
correct. Below we will demonstrate that on a perturbative level our resulta is equivalent to that of Ref. |^. 

Before coming to a more technical part of our analysis it would be useful to briefly remind the reader about the 
main steps of our calculation!! An attempt to examine our derivation "step by step" has been already made in 
Sec. 6 of Ref. |^. Unfortunately the information contained there is incomplete. In fact fromd the reader might get an 
impression that what we do essentially boils down to two trivial technical steps (performing the Hubbard-Stratonovich 
transformation and representing the Green functions in terms of the path integral) and one ("incorrect") disorder 
average. The actual procedurea is richer both formally and physically. 

Our analysisu consists of the two main steps. 

1. The first step is to reformulate the initial many-body problem with interaction in terms of a single quantum 
particle interacting with an effective quantum environment. Here we indeed use the Hubbard-Stratonovich 
transformation, but it is only a simple technical tool. A somewhat less trivial task is to derive a formally exact 
equation of motion for the single electron density matrix p in the presence of interaction (eqs. (24-25) of Ref. 

I): 

l ~W = ~ eV+lPv ^ ~ ( 1_ Pv)-^Pv - Pv^—(l - Pv), (Al) 



The matrix p is obtained after averaging of py over the fluctuating fields V± carried out with the (again formally 
exact) effective action S[V+, V_] derived by integrating out electronic degrees of freedom. No approximations 
have been made so far. 

Although the expression for S[V + , V-] is too complicated to deal with, some important observations can be made 
already at this stage. Namely, the fluctuating field V+ enters the equations just like an external field whereas 
the field V- enters in a qualitatively different manner. Fluctuations of the field V+ are essentially responsible 
for dephasing. 

In order to proceed further we make the first approximation: we evaluate the effective action iS[V+, V_] within 
RPA. The action S is now quadratic in V± and contains the dielectric susceptibility e(u>, k) of the effective 
environment. After that we easily integrate out the fields V± and arrive at the influence functional F for 
interacting electrons in a disordered metal. This completes the first part of our analysis. We can only add that 
this approach can be also used in the situations when approximations other than RPA are more appropriate. 
In those cases the action V_] and the influence functional should be modified accordingly. 

2. As a result of our derivation we arrived at the problem of a quantum particle in a random potential in the 
presence of the effective environment described by the influence functional F. The Fermi statistics and the 
Pauli principle are explicitely accounted for in the expression for F. The kinetic energy of a particle E is 
counted from the Fermi energy p and the states with E < are forbidden. The second stcEuof our analysis is 
to investigate the quantum dynamics of such a particle. According to the general principleaM at this stage the 
actual physical nature of the environment already plays no important role, any environment yields the same effect 
as long as it is described by the same influence functional. One can also develop better qualitative understanding 
of the phenomenon with the aid of simplier mojdels-.for the environment. In this sense the experience gained in 
the Caldeira-Leggett-type of models (see e.g.QB'Q~hj) is of particular interest and will be used in the Appendix 
B. Now we only mention one important feature of all these problems: the results for physically measurable 
equilibrium quantities do depend on the high frequency cutoff oj c of the effective environment. This dependence 
has nothing to do with the excitation of the environment oscillators, it exists even if the expectation values are 
calculated in the true ground state of the whole system. We will return to this point further below. 

Let us now come to a more technical part of this Appendix. We are going to test our expression for the influence 
functional at the level of the Golden-rule-type perturbation theory employed by the authorsa. Since at the first stage 
of our calculation we only used RPA (the same approximation was used inn) we have to recover all the same diagrams 
as in Ref. |[ However, AAG state that it is not the case. 
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We proceed in two steps. We first transform the AAG result for the conductance and demonstrate that by virt ue o f 
the c ausality principle one can completely remove the terms of the type G R G A G R G A . We will arrive at the Eqs. (|af| 
Af5) which are exactly equivalent to the resulta. Our second step is t o expand our expressio n for the conductance 
correctiond in the interaction terms. This will lead us to the Eq. ( A24 ) which is identical to ( A15 ). The reader not 
interested in technical details can skip the technical part and continue reading after the Eq. (A24). 

We startTrom reproducing the AAG expression for the correction to the conductivity due to electron-electron 
interactionO which they split into two terms &a a fi — ^ (T t^g h ^"afl > wnere 



c deph " 

*°J =-16 



dridr 2 dr^dr 4 
V 



dm de f d e 

— tanh — — 

2tt 2vr V de 2T 



coth — — 
2T 



tanh ■ 



2T 



{2j a [Gf 2 (e) - G A (e)} j p [G A (e)G A (e - u)G A (e) - G*(e)G* (e - u)G*{e)] 
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(A2) 



dridr 2 dr 3 dr 4 f dcu de ( d 



V 



duj de I d e \ 

— tanh — 

2tt 2tt V de 2T 



e — lu 

tanh x 

2T 



{j a [Gf 2 (e) - G A (e)] j fi [G 2 \(e)G A 1 (e) - G&(e)G£(e)] [G&(e - - G&(e - w)£&(u)] 

+ ... + a^(3} (A3) 

For simplicity we keep the same notations as ml: 

qR( a ) are the retarded (advanced) Green functions for noninteracting 
electrons and £ R ( A ^ are photon propagators. Here ... stands for the terms containing two w-dependent Green functions. 
Such terms do not contribute to 1/t v neither inEl nor in o ur analy sisa and therefore these terms will be ignored further 
below. For our purposes it will be convenient to rewrite ( A2, A3|) in the form 

dr 1 dr 2 dr 3 dr 4 f duj de ( d e \ w r rR , . r A , w 

~ J 2^ 27 [de tanh 2T J C ° th 2T ^ ] ~ C ^ * 



5<7 a R = — — 



16 



V 



{2j a [Gf 2 ( £ ) - G A (e)] j [G A (e)G A 4 (e - w )G&(e) - Gg(e)G*(e u)G*(e)] + ... + a «-> /?} 
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dr 4 dr 2 dr^dr 4 j" dw de ^ d ^ ^ < \ 



tanh 

2tt 2tt \de 2T 



tanh x 

2T 



{2j Q [Gf 2 (e) - Gf 2 (e)] & [G&(e)G£(e)Z&M - Gf 3 (e)Gf 1 (e)£f 4 (o;)] [G&(e - w) - G&(e - u,)] 

+... + a <->/?}. (A4) 

We observe that the factor tanh enters in this expresion together with the difference [G R 4 (e — u>) — G A 4 (e — u)] . 
This combination is just the Keldysh function 



G K (e,r 1) r 2 )=tanh^ ; [G R (e, n, r 2 ) - G A (e, r x , r 2 )] = tanh ^ 



1 



1 



e + // - H + iO e + fx - H - iO 



(A5) 



This function can also be rewritten as follows 
G K (e ) r 1 ,r 2 )=t&nh^=J2 



1 



e - £ A + «0 e - £ A - *0 



*A(ri)^(r 2 ) 



tanh £(-2™)<*(e - 6)*A(ri)*J(r 2 ) 



£(-2™) ( tanh eA)*A(ri)¥J(r a ) 



*A(n)^(r 2 ) 



= /" dr' [G i? (e,r 1 ,r')-G' 4 (6,r 1 ,r')] (<J(r' - r 2 ) - 2p(r',r 2 )), 



(A6) 



where £\> ^a are respectively the eigenvalues and the eigenfunctions of the Hamiltonian H — fj,; p(r',r 2 ) is the 
equilibrium single electron density matrix, p = l/(exp((H — p)/T) + 1). In a similar manner one obtains 
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(! tanh ^r) [G R (e,r 1 ,r 2 )-G A (e,r 1 ,r 2 )} =2 J dr' ^^p- [G R (e, r' , r 2 ) - G*(e, r', r 2 )] . (A7) 

We also introduce the evolution operator U(t) = exp(—i(H — fx)t) which is defined both for positive and negative 
times. The functions G R and G are related to this operator by means of the following equations: 



G R (t,r u r 2 ) = -i6(t)U(t,r u r 2 ); G A {t, n, r 2 ) = i0(-t)U{t, n, r 2 ). 



(A8) 



Now let u s wr ite down the two equivalent forms of the Keldysh Green function in the real time representation. Wc 
find from ( A5): 



+00 



G K (t,n,r 2 ) 



dt' 



-iT 



sinh(7rT(t - f )) 
T 



— [GV,n,r 2 )-GV,n,r 2 )] 



r/f' 



sinh(7rT(t - 1')) 



T^U{t', ri ,r 2 ), 



(A9) 



and from (A6) we get 



G K (t, ril r 2 ) = J dr' [G R (t, ri ,r') - G A (t, ri ,r')] (S(r' - r 3 ) - 2p(r',r 2 )) 

= -*/ dr'^(i,ri,r')(<J(r'-r 2 ) -2p(r',r 2 )). (A10) 

Analogously we obtain 

+00 

^tanh^) [ G H( e , ri , Pa )-^(e,r 1 ,r a) ] / ^-^^ [G R (t> , ri ,r 2 ) - G A (t> , ri ,r 2 )] 

-iT(t - 1') 



d£ 



sinh(7rT(t - t')) 



(All) 



and 



2 J dr' dp(r g ur,) [G R (e,r',r 2 ) - G A (e,r',r 2 )] ^2 J dr' 9 ^ 1 ^ [G R (t,r> ,r 2 ) - G A (t,r< ,r 2 ) 

,dp{r x ,r') 



= -2i / dr 



d/j, 



-U(t,r\r 2 ). 



(A12) 



It is easy to observe that the eqs. ( A^, All ) contain the integral over time which does not enter the eqs. ( A10 A12). It 
is this additional time integration that leads to violation of the normal time ordering at the level of the perturbation 
theory and is responsible for the appearence of the diagrams G R G A G R G A . The interpretation of such di agrams in 
terms of the path integral is not possible. However, if one uses the other form of the same functions ( |A10| , [A1^ ) 
the normal time ordering is automatically restored, the combinations G R G A G R G A dissapear due to the causality 
principle and the path integral interpretation of the remaining terms of the perturbation theory can be made. 

We emphasize that all the above transformations are exact and have the advantage that in the final expressions 
only the propagators depend on the frequencies e and uj (except for the factor coth ^ in 8a a p). This allows one to 
use the analytical properties of the propagators related to the causality principle. Namely, G R (e) and £ r (lu) have no 
singularities in the upper half-plane, while G (e) and C (iS) are analytic functions in the lower half-plane. Making 
use of these properties one can easily prove the identities 



(A13) 



duC R {u)G A {e - w) = 0, / deG A (e)G A (e)G A 4 (e - uj)G A (e) = 0, 
duj£ A (uj)G R (e - u) = 0, / deG R (e)G R {e)G R A (z - w)G£(c) = 0. 
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Consider e.g. the integral J duo C R {uj)G a (e — uj). Since both functions C R {uo) and G A (e — uj) are regular in the upper 
half-plane, the integral vanishes. Alternatively, we can write J du C R {uj)G a {e — uj) = 2n J dtexp(iet)C R (t)G A (t) and 
note that C R (t) = for t < due to the causality principle, whil e G (t) = for t > and the integral is identically 
equal to zero. Analogously one can prove all the other identities (A13). 
The corrections to the conductivity can now be considerably simplified: 



6a, 



a/3 



dr 1 dr 2 dr 3 dr 4: dr 5 f duo de uj { r a -, 



V 



Gf 5 (e)- 



dfi 



jpG A {e)G A A {e-")Gt{e)+j a 



Gf 5 (e)- 



jpG R (e)G R (e - uj)G R (e) + .. + a ^ 



dridr 2 dr 3 dr4dr^drQ f du de 
V / 2^2^ 



(e)^J<4(e) [G A (e - w) (1 - 2p) m ] G A (e)Cf 4 (uo) 



+ j a G A (e)^j G R (e) [G R (e - u>) (1 - 2p)J G R (e)£ R (uj) + ... + a <-> /?} 



(A14) 



We observe that the terms of the type j a G R 2 (e)j /3G 2 \(e)G 34 (e ~ uo)G A 1 {e) do not enter the expression ( A14 ) at all. 
For later purposes it will be useful to rewrite the above expression in the form of the time integral: 



6a, 



a/3 



where 
R(t,r) 

I(t,r) 



e 2 f dridr 2 dr 3 dr4dr$ 
~2 



V 



dt\ J dt 2 J dt 3 x 
o 







9p 5 



jaU^^'jpU+ih) [/34(*2 - t 3 )U+(t 2 - t 3 )] U+{tl - h) 

"jaUtSi)^h u 2i{h) [hiih - t 3 )U 34 (t 2 - t 3 )} f7 4 i(*i - ta) + ... + a 



le 



dr i dr 2 dr 3 dr 4 dr 5 dr $ 



V 



dti J dt 2 J dt 3 x 




-hUis{ti)^huUh) [Ru(t 2 - t 3 )U+(t 2 - t 3 ) (1 - 2 P ) 64 ] U+{t x - tt) 
+ j a U+{t l )^'jpU2z{h) [R 3i (t 2 - t 3 )U 36 {t 2 - t 3 ) (1 - 2p) 64 ] U 41 (h - t 2 ) + ... + a «-> 



dujd 3 k 47T 
(2tt) 4 fc 2 e(w,fc)' 



1 



£*(t,r) 



1 



dud 3 k 



(2*)' 



Tm 



-47T 



k 2 e(uo, k) 



coth f— ) 
\2T ) 



-iu)t-\-ihr 



2eH 



duod 3 k ( uj 



(A15) 



[£ fl (u;,fc)-/: A (c,fc)] 



Now we will demonstrate that the equation (A.15) can be obtained within the path integral formalism. The formal 
expression for the conductivity has the forma 



° = j ^ j dVildVi2 ( Vri / ~ Vr2 /) \r lf =r 2f J {~ t > t ''' r lfi r 2f'i r Ui r 2i)( r U ~ r 2i)Po(ru, T 2i ). 



(A16) 



The kernel J is given by the path integral over electron coordinates and momentums ri(f),pi(t) and r 2 (t),p 2 (t) 
corresponding respectively to the forward and backward parts of the Keldysh contour. The explicit expression for 
this kernel read 



ri(t)=r lf 



r 2 (t)=r 2 f 



J(t,t';rif,r 2 f,r 1 i,r 2 i 



Dri / Vr 2 J Dpi / Vp 2 x 

ri(t')=rii r 2 (t')=r 2i 

x exp {iS [ri,pi] ~iS [r 2 ,p 2 ] - iS R {r 1 ,p 1 ,r 2 ,p 2 ] - S I {r 1 ,r 2 }}; 



(A17) 
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where 



S [r,p] = J dt"(pr-^-U(r) 



(A18) 



e 2 f 

S R [r 1 ,p 1 ,r 2 ,p 2 ] = y / dti 


J dt 2 {R(h-t 2 ,r 1 (t l 


-ri(t 2 ))[l-2n(pi(t 2 ),r 1 (t 2 ))] 


V 

-R(h - 


t> 

h,r 2 (t x ) - r 2 {t 2 ))[l - 


2n(p 2 (t 2 ),r 2 (t 2 ))] 


+R(h - 


t 2 ,r 1 {h) - r 2 (t 2 ))[l - 


2n(p 2 (t 2 ),r 2 {t 2 ))] 


-R(h - 


t 2 ,r 2 {h) - n(t 2 ))[l - 


2n(p 1 (t 2 ),r 1 {t 2 ))]}; 



and 



(A19) 



Si[ri,n] = — J dti J dt 2 |I(ii-t2,ri(ii)-ri(*2)) + /(*i-t2,ra(t 1 )-ra(t2)) 

-J(ti - t a > ri(ti) - r 2 (t 2 )) - I(h - t 2 ,r 2 (h) - n(t 2 )) 
The functions R(t,r) and I(t,r) were already defined above. 



(A20) 



In order to obta in the perturbative result (A15) from the formally exact expression (A16) one needs to expand 
the kernel J ( |A17 ) in iSr + Si. In the first order one obtains eight different terms. Again we will consider only 
the terms contributing to 1/t v , i.e. the terms containing R(ti — t 2 ,ri(ti) — ri(t 2 )), R(t\ — t 2 ,r 2 (ti) — r 2 (t 2 )), 
I(t\ — t 2 ,ri(ti) — rx(t 2 )) and I(t\ — t 2 ,r 2 (ti) — r 2 (t 2 )). Diagrammatically, these contributions are described by the 
first two diagrams shown in Fig. 4 of Ref. 0. Four other terms which relate two different branches of the Keldysh 
contour and contain both n and r 2 , e.g. of the type I(ti — t 2 , ri(ti) — r 2 (t 2 )), are described by the last diagram in 
Fig. 4 of Ref. [| These ter ms gi ve the contributions containing two Green functions depending on the frequency w, 
and were denoted as ... in ( AI5| ). Such terms arc fully reproduced within our method as well, however we will not 
consider them here for the sake of simplicity. 

The correction to the kernel J due to the term I(ti — t 2 ,r%(ti) — T"i(t 2 )) has the form 



ri(t)=r lf 



r 2 (t)=r 2 f 



^"(MVy,^;^/,^/) = -e 2 J dt 3 J dt 2 J Vn J Vr 2 J Vp x J Vp 2 x 

V V r 1 (t')=r li r 2 (t')=r 2 i 

xl(t 3 -t 2 ,ri(t 3 ) -ri(t 2 ))exp{i5 [ri,pi] - iS [r 2 ,p 2 ]} 

t t s 

= — e 2 / dr-i / dr 4 / dt 3 / dt 2 x 



U r 2f ,r 2l (t - t')U rif , 3 (t - t 3 )hi{t3 - t 2 )U 3A {t 3 - t 2 )U 4 ,r u (t 2 - t'). (A2I) 

Here we made use of a simple property of a path integral: 

r(t)=r } 

I Vr J Vpf(t",r{t"))c^{iS [r,p]} = / dr"U(t - t";r f ,r")f(t" ,r")U(t" - t;r" (A22) 

r{t')—ri 



which holds for an arbitrary function f(t",r(t")). Actually in deriving (A2I) the property (|A22|) was used twice 
because the function of two arguments I(t 2 — t 3 ,n(t 2 ) — ri (t 3 )) enters under the integral ( |A21| ). Alrea dy at this 
stage one can observe the similarity between the expression (A2I) and the second term in the expression (A.15). To 
establish the equivalence between these two expressions the following steps are in order: i) after substituting the result 
( |A21j) into the expression for the conductivity ( |A16[ ) and applying the current operator j = (ie/m)(V ri , — V r2 ,) one 
puts rif = r 2 f = r 2 , ru = r\, r 2i = r 5 ; ii) one denotes t — t' — * t\, t — t 2 — * t 2 , t — t 3 — > t 3 ; hi) one introduces an 



G 



additional integration J dr 2 /V which is just averaging of the expression (A16) over the sample volume and iv) one 
transforms the effective initial density matrix as follows 

{rii - r 2 i)po{ri il r 2 i) = i 7 — —7 *Ai [ r ii)^\ 2 ( r 2z) - -1 t. ■ (A23) 

Ai A2 

After these transformations one can immediately observe the equivalence of the results obtained by means of two 
methodstl andEl on the level of the perturbation theory. The terms arising from the real part of the action Sr can be 
transformed analogously, the only differ ence in this case is the presence of an additional factor (1 — 2p)64 related to 



the term 1 — 2n(p, r) in the expression (A19). Finally we get 

+00 t\ t% 

S* = -jJ dr ^ dr ^ dr ^ J dtl J dt 2 J dt 3 x 





p ' '^uutSuUh) [hi(t 2 - h)U+(t 2 - i 3 )] U+(h - t 2 



m d/i 
P dp 15 
m djjL 

+ OO t\ t 

ie 3 f dridr 2 dr 3 dr 4 dr^drQ 



+ ^^U+(ti)3U 23 (t 3 ) [hi(t 2 - t 3 )U 34 (t 2 - t 3 )} U 41 (h - t 2 ) + .. 



dti J dt 2 J dt 3 x 
1: 



6 ./ V 

■U 5 2(h)jU+(t3) [#34(*2 - t 3 )U+(t 2 - t 3 ) (1 - 2p) 64 ] U+(h - t 2 ) 



P dPl5 ; 



P dpib 
m dfjL 



Ut 2 {ti)3U 23 {t 3 ) [R 34 {t 2 - t 3 )U 36 (t 2 - t 3 ) (1 - 2p) 64 ] U A1 {h - t 2 ) + ...}. (A24) 



In order to verify complete equ ivale nce of (|A15|) and (A24) one should a) replace the operator ep/m by j; b) adjust 



the factor 3 by observing that (A24) and flA15| ) are the corrections respectively to to the scalar conductivity and the 
conductivity tensor (in the i sotro pic case one has 6a = (5a xx + Sa yy + Sa zz )/3) and c) adjust another factor 2 having 



in mind symmetrization of ( A15p with respect to indices a and j3. Also, one should keep in mind that the ope rator 



dpis/ d/i commutes with the evolution operator L^. This completes the proof of equivalence of the results ( A15 ) and 



(A24) 



Thus the AAG statement that within our analysis "only some paths were selected by hand" is false. As it was 
explained above the paths presented e.g. in Fig. 10c of Ref. || cannot appear in the path integral, they are forbidden 
by the causality principle. Unfortunately AAG did not indicate the direction of the electron motion (e.g. by arrows) 
in their Fig. 10c. Otherwise it would be completely clear that at some parts of this path the electron moves backward 
in time. Such diagrams can appear from direct multiplication of the Keldysh matrices, but not in the path integral. 

In any case the technical issue with "missing diagrams" is settled and now we can come to the central question: if 
there are no calculational errors what causes the difference between the results of Refs. || and ||? As it was already 
discussed AAG proceed perturbatively in the interaction. They state that we do the same: according to AAG our 
procedure "is nothing but a perturbative expansion"!!! and our results are "purely perturbative"cl. If this were true, 
our final results could be compared directly indeed. However, this is not true. Of course, in some sense we also 
proceed perturbatively when we expand the exact effective action S^Vy in powers of V± in the exponent. But this is 
just RPA, and it has nothing to do with the Golden-rule-type perturbation theory developed irJj. In this sense our 
procedure is essentially nonperturbative and (within RPA) includes processes in all orders. 



On a slightly more formal level we can reformulate the difference as follows. The conductance ( A16) is determined by 



the path integral (A17). As it was demonstrated above the AAG perturbative procedure is equivalent to expanding 
J up to the first order in iSfj + Si, integrating over time in the infinite limits (this yields energy conservation) 
and averaging the result over disorder within the quasiclassical approximation. Our procedure does not involve the 
expansion in iSr + Si. We just evaluate the complete path integral within the quasiclassical approximation taking 
into account all saddle point paths and average over disorder. The latter average is quite trivial because the part of 
Sr which could be important for dephasing disappears already before averaging. 

One might think that both procedures should give the same result as long as the interaction is weak and the terms 
i Sr + Si are small. Here, however, one should be cautious: these terms are defined by nonlocal in time expressions 
( Al(j| , A20 ) and involve integrals over the two times. These terms are never small at low T as long as time is not 
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small. Therefore in general one cannot expand and should evaluate the whole path integral no matter how weak the 
interaction is. 

Thus not only the equivalence of the results of refs. and [| on the perturbative level but also a n im portant 
difference between these results is established. In hoth papers the conduc tanc e of a disordered metal (A16) in th 
presence of interaction was calculated. In our papera the full path integral (A17) was evaluated whereas the analysis 
is cajuivalent to the first order expansion of this path integral in iSa + Si. It is exactly due to this reason the results 
andtJ differ. The average over disorder emphasized intl is not important at this stage: the results are different already 
before averaging. 

We can also add that for simplicity in Ref. || we did not consider the effect of magnetic field H. This effect can be 
trivially incorporated into our analysis by adding the term erA(r) to the electron Lagrangian, A is the corresponding 
vector potential. The whole procedures remains the same except the action Sq wiU now depend on A. Averaging 
over disorder at the last stage of the calculation will yield the equations (2.42) fromH which define the dependenca-of 
the magnetoconductance curves on H. This part of the calculation is standard and was discussed in details e.g. ir£3. 
We would only li ke to emphasize that even in the presence of the magnetic field H one is not allowed to expand the 
path integral J ( A17 ) in powers of iSr + Si for all relevant fields th 3> r e , where th is the dephasing time due to 
the magnetic field, see eqs. (2.9) of Ref. ||. _ _ 

Having established the important formal difference between the proceduresO ancH we could, in principle, already 
conclude our consideration. At this stage one does not even need to introduce and discuss such concepts as "de- 
coherence" , "dephasing time" etc. Just viewing t v as a formal parameter extracted in a standard way from the 
magnetoconductance measurements one immediately obse rves that the perturbative calculation^ yields t v — > oo at 
T — > 0, while if one evaluates the whole path integral ( A17) (which contains all the perturbative termsEj plus infinitely 
many other terms which contribution is not smaller than that kept intl) one obtains a finite result at low T in agreement 
with many experiments. .We believe, however, that it is important not only to establish the formal difference between 
the two approachesQ andu but also to understand the physical reasons for this difference to occur. The corresponding 
discussion is presented below. 



APPENDIX B: 

Let us first recall the qualitative arguments of AAG concerning the effect of quantum decoherencei: "Consider 
a quantum particle which moves in the environment of harmonic oscillators: each oscillator is characterized by its 
frequency to. The result of the collision of the particle with an oscillator depends on the relation between hu> and the 
temperature." AAG first consider the case Tiuj <^T and find that in this case "the probability for the inelastic collision 
is substantial". Then they write: "The situation for hu> ^> T is quite different. Indeed, up to the exponentially small 
terms the oscillator is in the ground state and the particle has the energy smaller than hw" . Here we already have 
two questions: (i) it remains unclear whether the interacting or noninteracting ground state of the system is meant 
and (ii) it is also not clear whether the whole system is in equilibrium at all. But let us first finish the citation fromtl: 
"Therefore, the energy transfer is forbidden by the energy conservation and the collision is elastic. Therefore, there 
is no difference whatsoever with the collision with such an oscillator and with the quenched disorder, which definitely 
does not dephase" . 

Here we definitely agree with the last part of this sentence, namely quenched disorder indeed cannot dephase. The 
remaining part emphasizes the main physical difference betweend and our work. The AAG arguments apply if one 
considers a scattering problem in which the total energy of the system is a sum of those for the noninteracting particle 
and the oscillators and it (the energy) is conserved during the whole process. Then, if initially all the oscillators were 
in their noninteracting ground states and the particle energy was small, in the end this particle will have the same 
energy because none of the oscillators can be excited. Hence, the collision is elastic and no dephasing takes place. 

We, however, are interested in another physical situation. Namely, we would like to describe the properties of a 
quantum particle which is in equilibrium with an (infinite) bath of oscillators or close to this equilibrium. In this case 
the energy exchange between the particle and the bath is possible and the total energy of this interacting system is 
different (in our case larger) as comparedj-tQ susum of energies of noninteracting particles. Then for various physical 
situations one can demonstrate (see e.g.offirLLil) that interaction of a quantum particle with other quantum degrees 
of freedom is qualitatively different from the effect of a static potential. Since the system with harmonic oscillators is 
simple enough one can illustrate the difference between the two physical situations by means of an exact calculation. 
This will be done below. 

The density matrix of a quantum particle q interacting with other quantum degrees of freedom can be represented 
in the form 
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P(lf^'f) = J dq l dq / l J(q f ,q' f ,t;q u q' i7 0)p l (q l ,q' i ) } 
where Pi{qi, q[) is the initial density matrix and 



J = 



it 



Vqi I Vq 2 e^p{iS Q [qi(t)} - iS Q [q 2 (t)} - iS R [q 1 ,q 2 ] - Si[qi,q 2 ]). 



(Bl) 



(B2) 



Here So[q] is the action of a noninteracting particle, while. the terms S R /j account for interaction and define the 
influence functional for an effect ive quantum environmentd'El. It is easy to observe that the Eq. (B2) has essentially 
the same form as the Eq. (A17). If the environment consists of harmonic oscillators with frequencies u> n and unity 
masses, the influence functional can be easily found provided one assumes bilinear in coordinates interaction between 
the particle and the oscillators. Defining q + = (qi + q 2 )/2 and g_ = q\ — q 2 one readily obtainaZm 



C 2 /"* 

Sr = 22 — / dtidt 2 q^{ti)sh\{uj n (ti-t 2 ))q + {t 2 ), 



(B3) 



Si = ^T — c °th(^J J dt 1 dt 2 q-(ti)coa(uj n (t 1 -t2))q-(t2). (B4) 



where C is a constant which governs the strength of interaction. At this stage the problem (B2-B4) is qualitatively 
analogous to one we arrive in our analysis of interacting electrons in a disordered metal after we employ an RPA 
expansion of the exact effective action (cf. eqs. ( A17- A2(i| )). In a model with oscillators we deal with the quadratic 



effective action from the very beginning. 

Let us assume C to be very small and try to find the transition probability 



Wif = J dq f dq' f dq i dq' i i;f{qf)^}(q' f )J(qf,q' f ,t;q l ,qlO)ij l (q l )'ip*(q' l ) (B5) 

from some initial state wit h th e wave function tpi(q) an d the energy Ei to some other orthogonal state ipf(q), Ef just 
by expanding the kernel J (B2) in powers of C. The calculation is trivial (see e.g.H) and we only reproduce the result: 



W if =BeJ2 7rC 'y JdtJ dt' e -*(%-W) J (coth £ + l) [5(u n u ) 5(u> n + w )] (B6) 

Integration over t—t' in the infinite limits yields the delta-function 5(Ef+ui—Ei) which ensures the energy conservation. 
At T — > only positive u> contribute and the second delta-function in the square brackets fails. Thus in this case the 
transitions are only possible if Ef +lu u — Ei = 0, i.e. for Ei < Ef the transition never happens because the system has 
no energy to excite the oscillator u) n . Thus it always stays in its initial state, i.e. Wif = and the quantum coherence 
is never lost. At nonzero T and Ei < Ef we have Wif ^ 0, but only oscillators with small frequencies w n < T can 
take part in the transitions, so that Wif is small as long as T is sufficiently low. _ 

The above simple example just illustrates the qualitative arguments of the authpraj presented above. Exactly 
the same physical reasons are behind the cancellation of diagrams demonstrated intj. This cancellation will always 
take place at T — > within the Golden-rule-type perturbation theory just because of energy conservation. The only 
difference is that instead of the combination "coth+1" for Bose particles there appears the combination "coth — tanh" 
in the case of electrons. Thisjs well known and was also rederived from our formalism (see Sec. 5 of Ref. ||). 

Thus the AAG calculations is fully consistent with their qualitative arguments. Both describe the same physical 
situation: the total energy of the system "particle+oscillators" is equal to its value without interaction and fixed 
during the whole process. 

It is obvious, however, that the formalism allows to go beyond the simple Golden-rule-type perturbation theory and 
to provide a fu l l de scription of an interacting system. In order to demonstrate that let us perform a simple Gaussian 



integral in (B2-B4) exactly for the case of a free quantum particle with a mass m. Just for the sake of convenience 

of the oscillators 

[S(oj n - lu) - 5(uj n + lo)} = t]uj p , \w\<w c , (B7) 



we will assume a continuous spectrum of the oscillators 

irC 2 



E 



2u> n 







where oj c defines the high frequen^pfjaaJxiff and is a positive number. The exact results for the kernel J are well 
known and are presented elsewhereO'EruO. Here we only focus on the most important terms. One finds 



J oc exp 



- mfi(t)qti - mf 2 {t)(q-f - q- t ) 2 + 



(B8) 



where 7 = 77/771 and q±i/f are the initial/final values of q±, For the Ohmic bath ((3=1) one obtains 



did 



1 - e 



-2nTt 



h (t) = lJ*sJ*Jj ^coth-e-(-') =7T* + 7 ln. 27r(TK) 

-UJ a 



(B9) 



In the long time limit the function f 2 coincides with the equilibrium value of the average kinetic energy of the particle 

(E) EE (777472): 



h = (E) = 7 



dui lu coth 



2T 



2tt w 2 + 7 2 



7 dn^ 



2tt 



7 



T T 
— arctan — . 

7T 7 



(BIO) 



We observe that the particle q looses its coherence due to interaction with the bath of oscillators. Indeed in the long 
time limit we have fi(t) 3> f 2 and the kernel 



effectively reduces to 



(Bll) 



where L is the system size. Any perturbation of the density matrix will relax to the same equilibrium form 

p( qi ,q 2 ) = (l/L)e-^~^ 2 / L «, L 2 d = l/mf 2 (B12) 

which is not sensitive to the initial phase. As a res ult o f interaction with the bath even at T — > the off-di agonal 
elements of the equilibrium density density matrix ( |B12| ) decay on a typical length scale Ld ~ 1 / \Jr] ln(w c /7) set by 
interaction. The average value of the kinetic energy of the particle (E) ( |bTo| ) is not zero even at T — irrespectively 
to its initial energy. At sufficiently low T its value is also determined by the i nteracti on parameter 7. Another 
observation is that the high frequency cutoff uj c explicitely enters the expressions (BE.B1C) which in the limit T — > 
describe the true ground state properties of the system. The same results can be obtained within the imaginary time 
technique, or just by an exact diagonalization of the initial Hamiltonian of the system "particle+oscillators"Q. From 
the latter work it is particularly transparent that all the results depend on lu c just because all the oscillators (including 
the high frequency ones) "take part" in the diagonalization. By no means this implies excitation of such oscillators. 
Rather one can say that in the presence of interaction the noninteracting energy levels of the oscillators acquire a 
finite width and they can exchange energy with a particle in arbitrarily small portions. As a result of this exchange 
the particle energy is distributed as 



w(E) 



dp 
2?r' 



5 E 



V 

Ira 



f(p)<xexp(-E/2(E)), 



with the average value (E) given by eq. (BIO). Here we defined 



f(p) = L dq^p{q-)e' 



where p(q~ ) is the equilibrium density matrix ( B12| ) and g_ = q\ — q 2 . 

In order to avoid misunderstandings we would like to emphasize that we (on purpose) work in basis of "noninter- 
acting" eigenstates of the system. It is obvious that the full wave function of the total system as well as each of the 
eigenmodes obtained by an exact diagonalization always stay coherent. However, since the behavior of the particle 
q (and not that of the eigenmodes) is of interest for us, the reduced density matrix p(qi, q 2 ) should be studied. The 
decay of the off-diagonal elements ofpO on the length scale ~ Ld just implies that the bath in some conventional 
sense "measures" the particle positiontj. In principle the off-diagonal elements of p(qi,q 2 ) (and thus the coherence of 
the particle q) can be suppressed completely (Ld tends to zero if one e.g. chooses w c — > 00), while the eigenmodes 
of the total system obviously remain fully coherent. Thus it is quite useless to discuss the presence or absence of 
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quantum coherence in the interacting many-body system without discussing which quantity is actually calculated 
and/or measured in experiments. 

Let us also note that the p arti cular dependence of the results on the high frequency cutoff uj c is fully determined 



by the spectrum of the bath (B7). E.g. for f3 < 1 for most of the physical quantities of interest the cutoff does not 
enter at all, for /3 = 1 this dependence is logarithmic, while for f3 > 1 one gets a power-law dependence on u> c . In 
the latter case also the decay of perturbations of the density n matEix to its equilibrium value at T = is faster than 
logarithmic. All these features have been extensively studiedoBtTO and we will not go into more details here. 

We conclude that the presented above Golden rule approach fails to reproduce all the features obtained from the 
exact solution. E.g. at T = the Golden rule approach yields (E) = ryO = 0, while the correct result is given by eq. 
( B10| ) , a ccor ding to this approach the initial state does not decay as a result of interaction, while in reality it does 



(cf. (B8-B9)), according to (B6) no coherence can be lost at T = 0, while actually it is lost and even in equilibrium 
the off-diagonal elements of the density matrix decay on a typical length scale ~ Ld set by interaction. The Golden 
rule approach cannot correctly describe both the ground state properties of an interacting system as well as their 
low-lying excitations (cf. e.g. Ref. || where such properties were derived by means of an exact diagonalization). 

The example with a free damped quantum particle is not unique, of course. One can also consider e.g. a degenerate 
two-level system and observe that interaction with the CL bath leads to the effective renormalization of the transition 
amplitude A between the two levels^: A — ► A r = A(A/u> c ) T:r " , where a oc r\. In the weak interaction limit one finds 
A r = A(l — a ln(w c /A)). Again this is the property of the true ground state of an interacting system which cxplicitely 
depends on the high frequency cutoff lo c even for very small a. Again it is impossible to obtain the above result for 
A r within the Golden rule approach which at T = gives A r — A = aO = 0. 

Coming back to the problem of interacting electrons in a disordered metal we observe exactly the same situation. 
The Golden-rule-type perturbation theoryQ is not sufficient at low T just because it cannot correctly account for the 
ground state properties of an interacting system. Proceeding perturbatively one starts from a noninteracting system 
and imposes energy conservation. Obviously no dephasing can occur at T = because the electron energy (E) is 
zero in this case. But this simply means that one keeps the energy of an interacting system lower that it is in its 
true ground state. In the presence of interaction this energy is not zero (E) > even at T = 0. Although we do 
not calculate this energy explicitcly within our analysis it is obvious that the scale for (E) is set by interaction and 
the high frequency cutoff of the effective environment will enter. As in the above example of the exactly solvable CL 
model this dependence has nothing to do with excitation of high frequency oscillators. This is just the property of 
the true ground state of an interacting system. 

Thus we see that our problem has essentially the same qualitative features as the CL model. Of course, there exist 
also physical differences between these problems. For instance, there seems to be no quantity in the CL model which 
would be completely analogous to the decoherence time t v measured in the weak localization experiments. One can, 
however, qualitatively compare the length Ld in the CL model with the decoherence length L v ~ ^Dt v . Actually 
both quantities agree qualitatively if one establishes the correspondence between the interaction parameters and the 
bath spectra of these problems. More importantly, in both problems quantum decoherence appears as a result of 
energy exchange between the particle and the effective environment. [We emphasize again that we are working in the 
basis of noninteracting electrons which is obviously not the basis of the eigenstates in the presence of interaction.] In 
both problems the effect cannot be captured within the Golden rule approach where this energy exchange is forbidden 
at T = 0. 

In order to illustrate this point again let us "forbid" this exchange in our results by hand. Formally it implies that 
one should take the limit uj — > in the expression for the dielectric susceptibility of the effective environment in which 
case Im l/e(u>, k) vanishes. Obviously this procedure is not justified both p hysi cally and mathematically, but we just 



use it for illustration. For the CL model one can e.g. integrate over s — s' in (BE ) in the infinite limits, get 8{uj) and as 
a result recover only the first term jTt while the second will be missing. Then in the limit T — > one would observe 
no decay of the initial density matrix which would be incorrect. Analogously at T = one would get (E) = which 
would be incorrect too. The correct procedure is - as it was done above - to integrate over frequencies keeping the 
time finite and only then to send it to infinity. 

Analogously if one (just by hand) substitutes the time integral in eq. (71) of Ref. || by the delta function 8(uj) in 
order to provide the limit Im l/e(u>, k) — > (no energy exchange) one obtains 

1 e 2 r d d k 
T- 



t v a d J (27T) d fc 2 

The result is again zero at T = 0, but again the T-independent term wilLbe missing as in the above case of the CL 
model. [Note, however, that for large T one can recover the correct resulttJ proceeding in such a way. This is achieved 
if one cuts the integral at the lower limit at k ~ 1/L V ~ 1 / J L)tA . Although the above procedure is not rigorous it 
clearly illustrates the physical difference between the results^ andB. 
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APPENDIX C: 



In order to complete the critique of our theory AAG presented a collection of various experimental results and 
reported a strong disagreement with our theory. In some 3d systems the disagreement was found to be 4 to 5 orders 
of magnitude. This must look as a strong argument against our theory. 

However, looking closer at the experimental data one observes that in many systems chosen by AAG for comparison 
the mean free path is several or even one (!) Angstrom. Such strongly disordered systems are well beyond the 
applicability of our quasiclassical theory. Hence, for such systems the comparison cannot be carried out at all. On 
top of that, many systems have the granular structure. Our analysis cannot be directly applied to such systems as 
well, for instance because the simple formula for the noise spectral density 

(|K,, fc | 2 ) _ ^coth# 



7 3 — o 



a d k 2 



used in our theory works only for relatively small frequencies. [Here d is the effective dimension and a is a film 
thi ckne ss in 2d or a square root of a wire cross section area in Id.] In such systems the high frequency cutoff lo c 



in (CI) is determined not by the elastic scattering rate but rather by the effective capacitance of metallic grains. 
Examples will be analyzed in the first part of this Appendix. 

In the second part of this Appendix we will discuss several experiments in which the saturation of t v was observed. 
All these experiments are in a quantitative agreement with our theory. It appears that no alternative explanation of 
these experiments exists at the moment. 

It was always implied that our analysis is applicable for good metals with high conductivity and relatively low 



elastic scattering rate. Then and only then one can use a simple formula (CI) up to frequencies u> ~ l/r e . In such 
complicated and strongly disordered systems as granular percolating systems and metallic glass es o ne typically has 
l/r e ~ 10 15 -7- 10 16 Hz. This corresponds to energies ~ lO 5 !^ where the simple approximation (CI) does not work 
due to various reasons. In a general case the noise power spectrum is expressed by the formula 

which should be used for the comparison with experiments. Below^we will demonstrate by how much the noise can 
be reduced (and hence t v increased) in many systems discussed inO if one uses a realistic model for e{uj, k). For Id 
and 2d systems some possible sources of the noise reduction have been already analyzed in Sec. 4c of Ref. ||. 

Let us e.g. consider an experimental workEil where the decoherence time t v was measured in 3d Cu granular 
percolating systems. According totl the inverse inelastic time in this system is l/r e ~ 2 x 10 15 Hz and the elastic mean 
free path (extracted from the diffusion coefficient) is I ~ 4 x 10~ 8 cm. It is clear that our analysis cannot be applied 
to systems with such small values of I because the quasiclassical approximation should not really work there. But for 
a moment we ignore this fact a nd c oncentrate on the expression for e(u>, k). The authorsa use our formula (81) from 
Ref. U obtained from the Eq. (CI) and reported 5 orders of magnitude difference with the measured value for r c 



But the formula (CI) definitely cannot be applied for granular materials already at frequencies much lower than 10 15 
Hz! The effect of grain capacitances at not very low u> is crucially important and drastically changes the frequency 
dependence of e. Using the standard model for a granular metal one easily arrives at the conclusion that the high 



frequency cutoff lo c in (CI) never exceeds uj c ~ 1/RC and can be even smaller if the grain self-capacitance dominates 
over the intergrain capacitance (here C is the typical grain capacitance and R is the typical intergrain resistance) . In 
order to get a rough estimate for lo c in this case we make use of the typical grain size ~ (1 3) X 10 -6 cm and the 
resistivity p « 7 x 10~ 5 fl cm reported by the authorsEil. A reasonable estimate for R would be R ~ 100 Q. Although 
the grain capacitance is difficult to estimate from the grain size we can guess that the corresponding (renormalized) 
charging energy E c < Ec = e 2 /2C should not exceed E c ~ CL1 K (most probably Eq is even smaller because the 
effect of capacitance renormalization should be quite strong inEJ). This yields an upper limit estimate for the high 
frequency cutoff uj c ~ 10 11 -j- 10 12 Hz (most probably even lower). As a result the estimate for the decoherence rate 
1/t v extracted by AAG from our simple formula should be reduced by a factor ~ (w c r e ) 3 / 2 ~ 10~ 5 -j- 10~ 6 . We see 
that already such a simple estimate allows to completely remove 5 orders of magnitude disagreement reported by 
AAG for the experimento. An analogous conclusion can be reached in the case of other granular systems, e.g.E3. By 
saying that we do not want to claim a goo d agreement of our theory with these experiments. Rather we want to 



emphasize that a simple approximation (CI) (and hence our simple formulas derived on its basis) cannot be applied 
e.g. for granular systems at frequencies exceeding oj c <C l/r e . The whole theory should be substantially modified in 
this case. 

We believe that the same conclusion can be made for metallic glasses. There the problem to construct a reasonable 
model for the dielectric susceptibility appears to be even more complicated. Such systems are usually strongly 
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disordered and the elastic scattering rate is very high, typically l/r e ~ 10 15 4- 10 16 Hz. These values are of the same 
order as the conductivity and the plasma frequency. They are also comparable to the atomic interlevel spacing. It 
is not quite clear to which extent the glasses can be considered as homogeneous materials. In some cases granular 
structure was reported, so that some properties of metallic glasses can be similar to those of granular material s. T hese 
systems are clearly different from simple metals in many respects and therefore the simple approximation (CI) can 
hardly be applied already at moderately high frequencies, certainly well below 10 16 Hz. As in the case of granular 
metals the effective cutoff lo c ~ 10 12 4- 10 13 Hz would allow to completely remove the disagreement in all cases. But - 
we repeat - the systems are very complicated and it is not clear whether our formulas can be applied even at w < lu c . 

Another very important condition used in our theory is -ppl 3> 1 or, equivalently, tFT e ^ 1, where cf is the 
Fermi energy. Only provided this condition is satisfied the quasiclassical diffusion picture can be applied. Again in 
the experiments with 3d and 2d metallic glasses this condition is not satisfied. In all the experiments which cannot 
be described by our formulas the parameter €FT e is smaller than 5. Thus the simple quasiclassical approximation 
used in our paperuLdoes not work in these cases. An interesting illustration of this point can be extracted from the 
experimental dataHl. In this work the two 2DEG samples with €FT e ~ 3 and one sample with €FT e ~ 45 were studied. 
For the strongly disordered samples the measured t v at low temperatures exceeds the value predicted by our theory 
by the factor ~ 200-4 300, while for a weakly disordered sample the agreement between the theory and the experiment 
is very good (r£ x P(T min - 300mK) = 1.5 x lO" 10 sec, r* hoor (T = 0) = 2x 10~ 10 sec). 

Yet another problem is the accuracy of independent measurements of the parameters of such complicated systems 
as strongly disordered granular materials and metallic glasses. Many authors estimate the elastic scattering time r e 
using the value for the density of states obtained within the free electron model. Komori et al.Ea extracted the value 
r e for disordered Cu films directly from the measured magnetptjesistance curves and arrived at much higher values for 
r e than obtained by another method. If we use the estimatetS we will immediately conclude that our theory works 
well in the case of disoredered Cu films. Here we do not want to discuss which way of evaluation of r e is better. 
Rather we would like to emphasize that in many cases the accuracy of determination of the system parameters is 
clearly insufficient for quantitative comparison with theoretical predictions. . .. . 

It is also interesting to note that in the experiments with 3d metallic glassesESEZI the saturation of t v {T) was 
observed at relatively high temperatures, T ~ 1 4- 4 K depending on the experiment. According to the authors, 
this saturation can be explained neither by heating nor by the ef fe c t .of magnetic impurities. If so, what could be 
an alternative explanation for this effect? If, following the authoraijO, one adopts a free electron model of metallic 
glasses and estimates the values of er, r e and other parameters, one immediately arrives at the conclusion that the 
simple Drude formula (CI) is not applicable. in this case and the noise should be greatly reduced. If, furthermore, 
one assumes (again following the authorsEj'EZl) that in the metallic glasses the effect of electron-phonon interaction is 
more important than that of electron-electron interaction, then the experimental data should be compared with the 
eq. (82) of Ref. | The results are summarised in the Table I: 



Sample^-, l/r e , Hz p F l E F r e a, Hz c, km/sec l/r° xp (T = 0), Hz l/r°- ph (T = 0), Hz 

CuroAlsJbl 3 x 10 ib 4 2 2.9 x 10 ib 6 1 x 10 1,J 6 x 10 lu 

Cu 5 oY 5 «y 3.2 x 10 15 6 3 5.7 x 10 15 2 4 x 10 9 1 x 10 10 

Y 80 Si2(JlZl 1.1 x 10 16 2 1 1.8 x 10 15 2 1 x 10 9 3 x 10 11 

TABLE I. 

Here c is the sound velocity. Taking into account the uncertainity of the experimental parameters as well as the 
accuracy of the theory, we can conclude that the agreement between theory and experiment is good. 

The electron-phonon formulas also well describe the experiments in 2d Mg filmsE2l. The results are given in the 
Table II. 



Sample 


l/r e , Hz 


PfI 


Eft. 


(7, Hz 


c, km/sec 


1 /r cxp Hz 

/ 'tp,max; iizJ 


l/r°- ph (T = 0), Hz 


Mg3 


1.3 x 10 lb 


1.63 


0.8 


1.6 x 10 ib 


6 


9 x 10 u 


6 x 10 iu 


Mg4 


6.1 x 10 15 


3.55 


1.8 


3.6 x 10 15 


6 


7 x 10 9 


7 x 10 9 


Mg5 


2.16 x 10 15 


10 


5 


1 x 10 16 


6 


1.5 x 10 s 


5.8 x 10 8 



TABLE II. 



Again the agreement between the theory and the experiment is reasonable. Thus in the above cases our theory 
provides a natural explanation for the saturation of t v at low temperatures. No alternative explanation is known to 
us at the moment. 

The saturation of the decoherence time extracted from the magnetoconductance measurements was by now observed 
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in very many experiments. Typically it occurs already at relatively high temperatures of order 1 K. Although in some 
cases it can be due to extrinsic factors and/or magnetic impurities these reasons definitely cannot explain the saturation 
effect in all cases. Just for an illustration let us consider three experiments from those quoted in Ref. [|. 

The decoherence time saturation was observed iiE3. In Ref. ^ the authors argued that this saturation can be due 
to the effect of external noise. Let us not discuss this explanation but just accept it for a moment. AAG concluded 
that the external noise may effectively destroy the phase coherence without heat-Lag the sample if the resistance of 
the latter is much smaller than the quantum resistance 24 Kfi. In the expcrimentscil the resistance was smaller, so in 
that case the necessary condition was satisfied. I— . 

The decoherence time saturation was. also observed in the experiments^. In this case the resistance was at the MQ 
level, sq-that according to the authorsE3 external noise cannot be the reason for this saturation. At the same time the 
authors^ observed a rapid increase of the resistance which they interpreted as a crossover to Xke strongly localized 
regime. As this crossover occurs close to the temperature where r v starts saturating the authorsE3 concluded that this 
saturation is not meaningful because t v makes little sense in the strongly localized regime. Although this argument 
is essentially linked to the interpretation of the resistance increase as a Thouless crossover, and although even in this 
case we do not quite understand why the dephasing length L v extracted from the magnetoresistance measurements 
should stay constant in the strongly localized regime being 3 times smaller than the localization length, for a moment 
let us accept this argument of the authors too. I— . 

The third experiment we are going to discuss is one by Pooke et al£j. These authors also observed the decoherence 
time saturation in the temperature range T < 1 K in three samples with resistances 120, 240 and 360 KfL According 
to AAGl3 external noise cannot cause this saturation (unless it produces overheating) since the resistance of each of 
these 3 samples was much bigger than 24 Kf2. Although the authorsEj did not present the resistance curves one can 
hardly expect the Thouless crossover to take place in their measurements in the temperature range 0.1 -fl K: e.g. 
for the 120 Kf2 sample the dephasing length L v saturates at the level nearly 10 times smaller than the localization 
length. We believe that under these conditions it is already quite difficult to argue that in the experimentsEil the 
"saturated" part of the curve r v (T) is meaningless. Since the effect of magnetic impurities and heating ace. excluded 
in this case one should think about yet another explanation of the effect specifically for the experimentscJ. We are 
not aware of any proposal in this direction at the moment 

Mow, let us recall that the results of all 3 experimentsE3ErL-3 are in a quantitative agreement with our theory 
(seeac3). The parameters in these experiments differ by several orders of magnitude, but for practically all samples 
the agreement is within a numerical factor of order one. We believe this agreement can hardly be interpreted as a 
simple coincidence. Rather we can assume that this agreement strongly supports the validity of our calculation which 
offers the same explanation for these three and many other experiments and makes it unnecessary to search for a 
special explanation for each of them separately. In some cases the agreement is very good, for other experiments 
it is only qualitative, but this is rather a question of finding a proper model for e(w, k) and/or the high frequency 
cutoff tu c for the effective environment for this or that particular system. Different experiments are carried out on 
physically very different structures and it would be naive to expect that there can exist a unique formula which alone 
could quantitatively describe all existing weak localization experiments. But the saturation effect at low temperatures 
appears to be a unique intrinsic property of very many systems, and this effect is naturally explained by our theory. 

Finally, let us briefly discuss the experiments where, the crossover to the insulating behavior was observed. We 
believe that the experimental results in 2d structuresE3c3 quoted by the authorsU by no means can be used as an 
argument against our theory. In all these experiments the crossover to the insulating state was observed for strongly 
disordered systems (ei?r e was always of order one or only slightly larger). This is well beyond the applicability range 
of our theory. r~iri 

Furthermore, at least in some cases (e.g.E3'C3) the films could have the granular structure on the nanometer scale. 
If this is the case the insulating crossover can be also explained in terms of the Coulomb blockade of electrons on 
grains which does not reauire any phase coherence of these electrons at all. It is well known and was demonstrated 
experimentally (see e.g.c2l) that this effect occurs in granular metals and arrays at temperatures T < E c , where 
Eq is the effective (renormalized) charging energy of the grain. If the resistance of a granular array or film (or, 
equivalently, the intergrain resistance) is (roughly) above the quantum resistance unit R q ks 24 Kf2, the capacitance 
renormalization is not important and the bare charging energy Ec sets the temperature scale for the metal-to-insulator 
crossover T cr ~ Ec- The estimate Ec 10 K is resonable for the grain size in the nanometer range and is consistent 
with the experimental data. If the intergrain resistance is lower than R q the effective charging energy gets strongly 
renormalized due to charge fluctuations (see e.g.Eil). As a result one gets E c oc Ec exp(— AGR q ), where G is the film 
resistance and A is constant. For orderered square 2d granular arrays one has A = 2, and A is definitely bigger in a 
disordered case. This might also explain the exponential dependence of T cr on GR q reported e.g. in Ref. We can 
also quote recent theoretical results by Nazarovo who demonstrated that charging effects may persist in disordered 
conductors even without grains and tunnel barriess. All these facts as well as some other experimental details (like e.g. 
an obvious inconsistency of the experimental dataEa with the standard scenario of the orthogonal-to-unitary transition 
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in the magnetic field) suggest that interactions can play quite an important role in the above experiments. It is not 
our purpose to discuss this issue in more details here. But in any case it is quite clear that the above experimental 
data can hardly be interpreted as contradicting to our theory. _ . 

The same conclusion can be made concerning the experiments^ where the metal-to-insulator crossover in quasi- Id 
structures was reported and interpreted as the evidence for the Thouless crossover to the stronly localized regime 
at low T. It is clear, however, that the "noninteracting" scenario of this crossover cannot be applied in this case. 
For instance, at the crossover temperature T ~ 1 K the measured value of Tt v was found to be Tt v ~ 0.3 -j- 0.5 
depending on the sample and remained of order one even deep in the weak localization regime (e.g. at T ~ 10 K). 
As it was already discussed above the dephasing length L v was found to saturate at the value 3 times smaller than 
the "nonintecacting" localization length. All that implies-that interaction effects play a very important role in the 
experimentsEa. This was also acknowledged by the authors^. If so, the physical origin of the crossover can be debated. 
It is not clear how the Thouless scenario should be modified in the presence of strong interaction. There exist other 
mechanisms of localization which are solely due to interaction. E.g. at,T = a quantum particle can get localized 
in a periodic potential if it Interacts with an Ohmic bath of oscillators^!. In the subohmic case it is localized even in. 
the absence of any potentiate. Also the actual role of Coulomb-blockade-type of effects in disordered conductors^ 
should be understood better. All these effects do not involve the phase coherence of electrons and therefore would be 
compatible path the low temperature saturation of the decoherence length L v predicted by our theory and found in 
experiments^!. n 

Summarizing our discussion of the experiments we can conclude that in most cases discussed inEl the comparison with 
our theory cannot be carried out at all because the corresponding systems are very strongly disordered (the typical 
value of the mean free path is in the range I ~ 10 -8 cm) and therefore cannot be described by our quasiclassical 
theory. For other experiments the agreement is reasonable and can be improved further if a realistic model for e(u>, k) 
is chosen. In many cases, especially for quasi-ld systems, our theory is in a quantitative agreement with experiments. 
Further experiments are needed to study the effect of low temperature saturation of t v in more details. 
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